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Abstract 

We apply the quantum Singwi-Tosi-Land-Sjolander (QSTLS) theory for a study 
of many-body effects in the quasi-two-dimensional (Q2D) electron liquid (EL) in 
GaAs / Al^Gai-xAs heterojunctions. The effect of the layer thickness is included through a 
variational approach. We have calculated the density, spin-density static structure factors, 
spin-dependent pair distribution functions (PDF) and compared our results with those of 
two-dimensional (2D) EL given in earlier papers. Using the static structure factors (SSF) we 
have calculated various dynamic correlation functions such as spin-dependent local-field factors 
(LFF) and effective potentials of the Q2D EL. We have also calculated the inverse static dielectric 
function of the 2D and Q2D EL using different approximations. We find that the effect of finite 
thickness on the dielectric function is remarkable and at the intermediate values of wave number 
q there is a significant difference between the QSTLS and STLS results. 

PACS numbers: 73.21.Fg, 71.45. Gm, 71.10.-w 
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I. INTRODUCTION 



During the last decades quasi-two-dimensional electron systems have been of consider- 
able interest because of technological relevance to high-mobility electronic devices. Many 
authors have investigated Q2D systems applying various methodsi and most of the theoret- 
ical calculations have been performed in the framework of the random phase approximation 
(RPA). However, it is well-known that the RPA neglects the Pauli and Coulomb hole sur- 
rounding each particle and is not satisfactory approximation at low densities^. Therefore 
a number of authors have studied beyond-RPA effects by including the static LFF such 
as in the Singwi-Tosi-Land-Sjolander (STLS) approach'^. A common viewpoint is the need 
to incorporate dynamic correlations and the importance of dynamic LFF is evidence. The 
dynamic LFF has been introduced in several works via different schemes^. One of these 
works is that of Hasegawa and Shimizu using the Wigner distribution function-. The ap- 
proach is usually named as the quantum STLS and has been applied to 2D EL by Moudgil 
and coworkers^. Their calculations have been generalized in our recent work^ to the more 
realistic case of electrons in Q2D semiconducting heterostructures investigated by Bulutay 
and Tomak (BT)^. The aim of this work is to extend our previous investigations to include 
both electron and spin correlations. 

The outline of this paper is as follows. In Sec. II, we discuss briefly the theoretical for- 
malism of the self-consistent QSTLS equations applied to the Q2D EL in GaAs/Al^.Gai_2;As 
hetero junctions based on the variational approach proposed by Bastard^. In Sec. Ill the 
results and discussion are given, and in Sec. IV we summarize the main results and gather 
our assessment of the performance of the QSTLS for semiconductor heterostructures. 

II. THEORY 

We consider a GaAs/Ala;Gai_a;As heteroj unction where the electrons move into the GaAs 
side and form 2D subbands. For an accurate account of the electronic distribution we use 
Bastard's variational approach given in BT's work, where the electron effective masses mA,B 
and dielectric constants eA,B are considered to be different in the GaAs and Ala;Gai_a;As 
layers. In the regime where only the lowest subband is populated the self-consistent equations 
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for density and spin-density static structure factors at zero-temperature can be written as^ 
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with /° {p) is the Fermi - Dirac distribution function, xo [q] ^) = Xo (?, q] ^) is the free- 
electron response function and 

V (g) = (g) (5) 

qt 

is the Fourier transform of the effective potential. Here F{q) is the form factor to the 
Coulomb interaction due to the layer thickness given in the BT's work and £= {ea+Sb)/'^ 
is the average dielectric constant. A strictly 2D electron gas with (5-function density dis- 
tribution is obtained by setting F{q) = 1 in Eq. (5). To solve the sets of Eqs. (la-3a) 
and (lb-3b) we have used the procedure proposed by de Freitas et alM and obtained the 
following self-consistent equations 
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Here n, kp = \j1tui and are the electron density, Fermi wave number and angle between 
g and g*, respectively^. Using the integration by parts we can write the density and spin- 
density static structure factors as 
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These forms of the structure factors are appropriate for the numerical integration because 
the integrands remain finite in the limit (3 —>■ 0. 



III. RESULTS AND DISCUSSIONS 

In our recent papei^ we have solved the equations (6a- 7a) for 2D and Q2D electron liquids 
in semiconductor heteroj unctions. In the case of 2D EL we have obtained the SSF and PDF 
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in good agreement with Monte Carlo results. In this paper we solve the equations (6b- 7b) for 
the Q2D electron liquid in GaAs/Al2;Gai„^As heteroj unctions having a step-barrier potential 
with Uh =0.3eV , Ea =13, Eb = 12.1, =0.07me and =0.088me , where rrie is the 
vacuum mass of the electron^. Using the obtained results we calculate the spin-dependent 
pair-correlation functions, dynamic local-field factors, spin-dependent effective potentials, 
inverse static dielectric function and compare our results with those given in Refs.7 and 8. 



A. Static spin-density structure factors 

The static structure factor S{q) of the 2D EL was shown in Fig. 1 of our recent paper—. 
In this work we have calculated the spin-density SSF S{q) of 2D and Q2D EL by solving 
the Eqs. (6b) and (7b) in the self-consistent way for several values of electron density and 
the results are shown in Fig. 1. It is seen that in the case of 2D EL our results are similar 
to those of Ref. 6 and the effect of the layer thickness is remarkable for a wide range of 
electron densities. 



B. Spin-dependent electron pair-correlation functions 

The spin-symmetric and spin-antisymmetric electron pair-correlation functions can be 
expressed as 

(7n(r) = 0.5 [g{r)+g{r)] (13) 
gu{r) = 0.5 [g{r)-g{r)] (14) 

where 
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oo 
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To study the effect of electron density we show in Fig. 2. the spin-dependent PDF gi^{r) 
of 2D EL for different values of density parameter r^. We observe that our PDFs differ from 
those of Moudgil and coworkers at small values of r. This difference in the behavior of PDFs 
stems from the incorrect results of Moudgil and coworkers for the SSF at large q's discussed 
in our previous paper—. 
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To study the effect of the layer thickness on PDFs we have calculated the spin-dependent 
electron pair-correlation functions gii{r) of 2D and Q2D EL for different values of electron 
density parameter r^. We find that gii{r) is almost independent of and is therefore plotted 
in Fig. 3 only for Vg = 3. We observe from the figure that the effect of the layer thickness is 
considerable for an intermediate region of the inter-particle distance. 

C. Spin-symmetric and spin-antisymmetric dynamic local-field factors: 

To compare our results with those of Moudgil et al. and to study the effect of the layer 
thickness we show in Figs. 4 and 5 the spin-dependent dynamic LFFs of 2D and Q2D EL 
as a function of a; for q = l.lkp and = 3. We observe that our results for 2D EL are 
similar to those given in Ref.6 and the layer thickness has a remarkable influence on G(q,uj) 
and J(q,uj) for a wide range of frequencies u). 

D. Effective dynamic potentials 

Spin-symmetric and spin-antisymmetric dynamic effective potentials can be calculated 
from the spin-symmetric and spin-antisymmetric dynamic LFF as 



The effective potentials of the 2D and Q2D EL obtained in the static limit ( i.e., = 0) 
for = 1 and 3 are shown Figs. 6 and 7. We observe remarkable differences in the results 
of 2D and Q2D EL for all values of q and r^. We note that because of incorrect results for 
the SSF given in Ref. 6 our results for both spin-symmetric and spin- antisymmetric static 
effective potentials of 2D EL differ considerably from those of Moudgil and coworkers in the 
region of large values of q {q yS.bqp)- 

Following the authors of Ref.6 we define the spin-dependent effective dynamic potential 

as 
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The real and imaginary spin-dependent effective dynamic potentials of 2D and Q2D EL 
for = 1 and q =l.lfci? are plotted, respectively, in Figs. 8 and 9. It is seen from the 
figures that the effect of layer thickness is significant for a wide range of frequency uj and our 
results for 2D EL are similar to those of Moudgil and co-workers^. We note that the values 
of the real spin-dependent effective dynamic potentials shown in the Fig. 9(a) of Ref.6 is 
not correct. By using the values of V/^j {q,uj) and V^^j (q,^) at = we can see that the 
authors of Ref.6 have divided the correct values of cj) and VJ=//||(g, cj) by 2. 

E. Inverse static dielectric function 

Finally, we calculate the inverse static dielectric function of 2D and Q2D EL for =3 
using different approximations. The figure 10 shows that the QSTLS results differ less 
from those of the STLS approximation in the Q2D EL than in the 2D EL. However, at 
the intermediate values of q there is a significant difference between the QSTLS and STLS 
results. 

IV. CONCLUSIONS 

Using the QSTLS approximation we have studied the electron and spin correlations in 
GaAs/A^Gai_2:As heterojunctions. We have calculated the spin-dependent PDF, SSF, 
dynamic LFF, effective potential and inverse static dielectric function of 2D and Q2D EL. We 
have shown that the results for Q2D EL in semiconductor heterostructures differ remarkably 
from those of 2D EL and the difference between the QSTLS and the STLS results of Q2D 
EL is considerable. It is hoped that our results will be of help in investigating the effect of 
electron and spin correlations on properties of Q2D electron systems. 
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FIG. 4: Spin-symmetric dynamic structure factors for rg = 3 and q =l.lkF- 
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FIG. 9: Imaginary spin-dependent effective dynamic potentials = 1 and q =l.lkF. 




FIG. 10: Inverse static dielectric function of 2D and Q2D EL for =3. 
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